Absence of fermionic quasi-particles in the superfluid state of the attractive Fermi gas 
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We calculate the effect of order parameter fluctuations on the fermionic single-particle excitations 
in the superfluid state of neutral fermions interacting with short range attractive forces. We show 
that in dimensions D < 3 the singular effective interaction between the fermions mediated by 
the gapless Bogoliubov-Anderson mode prohibits the existence of well-defined quasi-particles. We 
explicitly calculate the single-particle spectral function in the BEG regime in D = 3 and show that 
in this case the quasi-particle residue and the density of states are logarithmically suppressed. 
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Recent experimental progress in cooling atomic Fermi 
gases to ultracold temperatures [I| has revived the in- 
terest in the BCS-BEC crossover in the attractive two- 
component Fermi gas. Although a simple mean-field ap- 
proximation is sufficient to obtain a qualitatively cor- 
rect description of the thermodynamics of this crossover 
@i S| I for comparison with experiments numerically ac- 
curate calculations beyond the mean-field approximation 
are necessary. Fluctuation corrections to the free energy 
have been calculated at the level of the ladder approxi- 
mation which in the functional integral approach is 
equivalent to the Gaussian approximation in the order 
parameter fluctuations [S i, 0, i, i, E El H H [11 [H ■ 
Close to the unitary point, where the s-wave scattering 
length is large, the Gaussian approximation is not suf- 
ficient and more sophisticated many-body methods are 
necessary [lil [l7t. 

While most authors focused on the thermodynamics of 
the BCS-BEC crossover, not much attention has been 
payed to the fermionic single-particle excitations in the 
BEC regime. Fluctuation corrections to the fermionic 
self-ene rgy have been briefly considered in the appendix 
of Ref. [iSj , but the momentum and frequency dependent 
single-particle spectral function p{k,Lo) and the density 
of states have not been analyzed. In this Letter we 
shall show that in dimensions D < i the singular inter- 
action between neutral fermions mediated by the gapless 
Bogoliubov-Anderson (BA) mode leads to the breakdown 
of the quasi-particle picture. We emphasize that our 
result is relevant for superfluid gases of cold fermionic 
atoms, but does not apply to the superconducting state 
of charged fermions, where the long-range Coulomb in- 
teraction pushes the B A mode up to high energies • 

For our explicit calculations of p{k,u)) and v{ljj) we 
shall focus on the BEC regime in the physically most 
relevant dimension D — i aX zero temperature. It turns 
out that in this case the imaginary part of the self-energy 
vanishes linearly in the properly shifted frequency vari- 
able y (see below), implying that the frequency depen- 
dent real part of the self-energy vanishes as ylny. A 
similar frequency dependence of the self-energy in the 



normal state of strongly correlated electrons has been 
proposed phenomcnologically by Varma et al. jlij to ex- 
plain the unusual normal state properties of the high- 
temperature superconductors. While a generally ac- 
cepted theory explaining the microscopic origin of such a 
marginal Fermi liquid behavior in the normal state does 
not exist, we shall give here a microscopic derivation of 
similar marginal quasi-particle behavior of the fermionic 
single-particle excitations in the superfluid state of neu- 
tral fermions in three dimensions. 

To begin with, let us briefly outline the derivation of 
fluctuation corrections to the fermionic self-energy using 
fimctional integration. We consider a system of neutral 
fermions with dispersion — k"^ /{2m) and a short-range 
attractive two-body interaction gp > depending on the 
total two-particle momentum p. Writing the fermionic 
single-particle Green function G in the Nambu-Gorkov 
basis [18| as a functional integral and decoupling the in- 
teraction in the spin-singlet particle-particle channel by 
means of a complex bosonic Hubbard-Stratonovich (HS) 
held '0J we can represent G as a functional average. 



G 



/P[(5V',5V^]e-^''«[*'^''^>lG^ 



(1) 



where G^ is the fermionic Green function for a fixed 
configuration of the HS field. Its inverse is given by 
the following matrix in Nambu-Gorkov and momentum- 
frequency space [20|, 



-IpK'^K Sk.K' (it^ + C-fc) 



(2) 



Here £,k = ~ f- where is the chemical potential, and 
the effective action S'^ff [(5-0, 5-0] of the HS field is 



glr[^^ (3) 

n=2 



where 6^p = ipp — ^pqAq, and the mean-field Green 
function Gq is obtained from Eq. ^ by approximating 
the HS field by its saddle point V-'p ~ Sp^Ao. The fluctu- 
ation matrix V is defined byV = Gq"^ — Gl^. The saddle 
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point condition leads to the usual BCS gap equation for 
Aq and the quasi-particle dispersion Ek — + Aq. 
For convenience we neglect the momentum dependence 
of 3p ~ Sq regularize the resulting ultraviolet di- 
vergence [7| in the BCS gap equation by eliminating go 
in favor of the corresponding two-body T-matrix in vac- 
uum, which we denote by g. The relevant dimensionless 
interaction is then g = v^g, where vq is the density of 
states at the Fermi energy in the absence of interactions. 
In Z? = 3 we have g — —2kFas/TT, where as is the s-wave 
scattering length in vacuum and kp = mvp = ^/2rnep is 
the Fermi momentum. 

The exact fermionic propagator G defined in Eq. ([1]) 
can be written as G = [Gq^ — where the self- 

energy matrix S involves a normal component '^{K) and 
an anomalous component SA{K) (which can be viewed 
as a fluctuation correction to Aq), 



\KK' 



Y.{K) SA{K) 
SA*{K) 



(4) 



The leading contribution to S can be written as S = 
(VGqV), where (...) denotes functional average with 
the effective action ([3|) in Gaussian approximation, re- 
taining only quadratic terms in the fluctuations. As 



emphasized by Castellani et al. [2l| (see also Ref. |22l|). 
the scaling behavior of the order parameter correlation 
functions is more transparent if we express the com- 
plex field Sipp in terms of two real fields [i^ by setting 
^'>pP = [xp + i4>p]/V^ and iJi/ip = [x_p — i0_p]/\/2. The 
fields xp and (j)p correspond to longitudinal and trans- 
verse fiuctuations of the superconducting order param- 
eter, respectively. The effective action ^ in Gaussian 
approximation is then [20| 



Ses[x,4'] ~ \jjyX*p,' 



r p r p 
r p r p 



XP 

4>p 



, (5) 



where for small p and lo the matrix elements of the 
bosonic propagator are 



pX0 

^ p 



^ p 



^ P 



i2C P 



Dc^ 



(2Ao)2 



^2p2 ■ 



2Aou) 



f,2p2 



Introducing the dimensionless integrals 



Ao^fc 



(6a) 
(6b) 
(6c) 

(7) 



we have Zi = hh/ilf + /|] and Z2 = [I^ + 

while the velocity c of the BA mode can be written as 

Dc^/vp = Ji/[/f -|-/|]. Numerical results for Zi, Z2 and 




FIG. 1: (Color online) Factors Zi, Z2, and c/vp defined in 
the text as a function of —1/g — n / (2kFas) for D = 3. 



c/vp in D = 3 as a function of the relevant dimension- 
less coupling —g~^ are shown in Fig. [T] In _D dimen- 
sions we obtain in the BCS limit c ~ vp/^/D, and in 
the BEC limit c — \J PBU^ImB, with the effective in- 
teraction t/* — ^^\g\- Here pB = p/2 and itib = 2m 
are the density and the mass of bosons formed by paired 
fcrmions with density p and mass m. Our result for c in 
the BEC limit agrees with the well-known Hartree-Fock 
result for the velocity of the elementary excitations in the 
interacting Bose gas provided we identify C/* with the ef- 
fective Hartree-Fock potential. Note that only in three 
dimensions C/* is given by the zero-energy T-matrix |(?| of 
the underlying fermionic system. In the limit Z? — > 2 the 
T-matrix vanishes as g ^ ~{D — 2)/vq (see Ref. |2j|), so 
that in two dimensions J7* = 2/i'o and hence c = vp/y/2 
for all values of go [2^ . Carrying out the Gaussian aver- 
age in S = (VGoV), we finally obtain for the fermionic 
self-energies defined in Eq. ([4]), 



S(i^) = -^l^[Ff + F^^ + 2tFf]Bo{P-K), 



(5 A (AO 



[Ft.^-F^^]AoiP-K), 



(8) 
(9) 



where Bo{K) = -[ico + ^k]/[uJ^ + and AoiK) = 
— Ao/[cl'^ -I- E^] are the normal and anomalous compo- 
nents of the mean-field propagator Gq. 

It turns out that in the BEC limit, where the dimen- 
sionless parameter A = 4mc^/Ao « Ao/{2\p,\) is small 
compared with unity (see Fig. [2]), the integrations in 
Eqs. (|8l9p can be explicitly carried out in T) = 3 at zero 
temperature, provided the effective interactions mediated 
by the bosonic fields are approximated by Eqs. ((6aH6c)) . 
In this regime the s-wave bound state energy « 2\fi\ is 
the largest energy scale in the problem. Then we may ex- 
pand the mean-field dispersion as E^ ~ i?o + fc^/(2TO»), 
withi;o « \fj,\ + Al/es = \p\{l + 2\'^) andm/m, w l-2\'^ 
for small A. In fact, the leading non-analytic behavior of 
the self-energies is entirely due to the infrared singular- 
ity in the propagator Fp'^ of the BA mode in Eq. (fFb]) . 
so that we may neglect the contributions from Fp'^ and 
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FIG. 2: (Color online) Dimensionless parameter A — 
4mc^/Ao and A — 8\^/tv in D = 3 as a function of —l/g. 
For comparison we also show the ratio Ao/(2|/ij), which ap- 
proaches 4mc^/Ao in the BEC limit. 



Fp^ in Eqs. (|8|9p . Within this approximation, which is 
accurate for \k\as ^ 1 and \ \uj\ — E-k\ <C we obtain 

Y:(k,iuj) = -2Ao[X^C{x,z)-{l + X'^)C{x,z*)], (10) 
5A{k,iuj) = 2AaX[C{x,z) + C{x,z*)], (11) 

where x = \k\/kc and z = (iuj — Ek)/ujc^ with kc = mc 
and ll>c = m(? 12. The complex function C{x, z) is 



TT C{x, z) 
■1 



X In 



2 + ln(-)--{y(TWT^ 

x)\ (12) 



\ + x- ^(l + a;)2-|-z 



(x 



where the ultraviolet cutoff M « (kcO-s) ^ takes into 
account that in Eqs. ()10|lip we have used the long- 
wavelength approximation (j6bp for . For small x and 
|z| the asymptotic expansion of C(x^z) is 

7rC(x,z) ^ ln(AfV4) + (^/2)ln(-l/z)-hO(z,.x2),(l3) 

so that after analytic continuation z ^ y + iO we obtain 
for the imaginary part ImC(x,y -I- iO) ~ Q{y)y/2 for 
small X and y. The non-analytic term z In 2; in Eq. (jl3p 
can also be derived directly from Eqs. (|8l9p by simple 
power- counting, which reveals for £) < 3 an even stronger 
algebraic singularity. Although our explicit calculation is 
only valid in the BEC regime, the power counting analy- 
sis of Eqs. (|8l9p shows that the coupling of the fermions 
to the gapless BA mode prohibits the existence of well- 
defined quasi-particles for Z) < 3 in the entire range of 
the BCS-BEC crossover. This nicely fits to the fact that 
in the interacting Bose gas the Bogoliubov fixed point is 
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unstable in dimensions -D < 3, see Ref. 

Given Eqs. (jlOllip . we may calculate the nor- 
mal component of the spectral function p{k,uj) = 
—TT^^lmB{k,uj + iO), where B{K) = B{k,iuj) is the up- 
per diagonal element of the matrix-propagator G defined 
in Eq. From the derivation of Eqs. (I10|ll|) it is clear 
that these expressions are only valid for |fc| <C •\/2r7i|/i| 
and ||w| — <C 1^1 . However, due to the hierarchy of 



energy scales Wc ^ Aq <C |a<| in the BEC regime, there 
are three characteristic regimes where the spectral line- 
shape exhibits rather distinct behavior: 

Regime I: Aq ^ maxje^, ||ti;| — Ek\} ^ In this 
regime the spectral function exhibits well-defined quasi- 
particle peaks at a; = ±i?fc, both of which have approx- 
imately the same height. The quasi-particle damping is 
7fc = 2esA'^(|fe|as)^, which is small compared with Ek, 
so that the fermionic single-particle excitations are well- 
defined quasi-particles. 

Regime II: ujc maxjefc, | |w| — <C Aq. In this 

intermediate regime the weight of the negative frequency 
peak is a factor of A^ smaller than the weight of the 
positive frequency peak, which agrees with the Hartree- 
Fock result. There is a considerable asymmetry between 
the positive and the negative frequency part. The quasi- 
particle damping in this regime is small and momentum- 
independent, 7fc w 2AoA^. 

Regime III: maxjefc, ||a;| — Sfc]} <C uic- Here it is nat- 
ural to measure momenta and energies in units of the 
natural scales kc and ujc associated with the BA mode. 
Writing p{k,uj) = 9(cj)p+ -I- 6(— u;)p„, the positive and 
negative frequency part of the spectral function can be 
written in scaling form, p±{k,uj) = ujc^p±{x,y), with 
X = |fc|/fcc and y = (|w| — Ek)/ujc. Graphs of the scal- 
ing functions p±{x,y) are shown in Fig. [3l Note that in 
regime HI, where x and \y\ are both small compared with 
unity, there is no spectral weight for \uj\ < Ek. More gen- 
erally, for negative y and x < 1 + \/\y\ the spectral func- 
tion vanishes because Im C{x, y + iQ) = in this regime. 
Physically, this is due to the kinematic constraint that a 
quasi-particle with momentum k can only spontaneously 
decay into another quasi-particle and a BA quantum if its 
velocity |fe|/m* exceeds the velocity c of the BA mode. 
The scaling functions for x,\y\ <C 1 can be written as 
P+{x,y) = p{y) and p-{x,y) = [A^ + {y/%Y]p{y), where 

p{y) = Q{y)iT-^Z{y)-i{y)l[y^ + -f^y)], (14) 

with Z{y) = 1/[1 — Alny] and the damping function 
^{y) = 7rA?//[l — Alny]. The effective coupling parame- 
ter A = BA'^/tt is shown in Fig. [5] as a function of — 
For Lo Ek the damping 7(y) is only logarithmically 
smaller than the real part of the quasi-particle energy 
while the quasi-particle residue Z{y) vanishes logarithmi- 
cally. Such a behavior resembles the marginal Fermi liq- 
uid scenario postulated phenomenologically in Ref. [l9j . 

The spectral function implies a logarithmic re- 
duction of the density of states ^{ijj) = J -^^p{k,Lu) for 

Lo Eq w I/^I + Aq/es. For positive frequencies and 
X\n[uJc/{oJ — Eq)] ^ 1 we obtain 



1^0 



A In 



LU - En 



(15) 



The first factor is the usual square-root singularity in 
D = 3. For uj — > —Eq, the asymptotic limit of i'{uj) is 
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FIG. 3: (Color online) Scaling functions p±{x, y) of the single- 
particle spectral function in the regimes II and III defined in 
the text. The plots are for A = 8lJc/Ao =0.05. In the regime 
III both X and \y\ are small compared with unity. For better 
comparison with p+{x,y) we have divided p-{x,y) by A^. 

a factor of smaller than for positive frequencies and 
involves the same logarithmic suppression. 

In summary, we have shown that in I? < 3 the coupling 
of neutral fermions to the gapless BA mode prohibits 
the existence of well-defined fermionic quasi-particles in 
the superfluid state. Although our explicit calculation is 
only controlled in the BEC limit, we have argued that 
it remains qualitatively valid for the entire BCS-BEC 
crossover, because the non-analyticities leading to the 
breakdown of the quasi-particle picture are due to phase- 
space restrictions in the infrared regime. In particular, 
close to the unitary point 1/g = 0, where all energy scales 
are of the order of , we predict a logarithmic suppres- 
sion of the density of states at the scale ep which should 
be observable via tunneling experiments or other exper- 
imental probes of the density of states in ultracold gases 
of neutral fermions. 

This work was supported by the DFG via 
SFB/TRR49. 



[2] D. M. Eagles, Phys. Rev. 186, 456 (1969). 

[3] A. J. Leggett, in Modem Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski and R. Przystawa, 
Lecture Notes in Physics 115 (Springer, Berlin, 1980). 

[4] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 
195 (1985). 

[5] M. Drechsler and W. Zwerger, Ann. Phys. (Leipzig) 1, 
15 (1992). 

[6] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 
Phys. Rev. Lett. 71, 3202 (1993). 

[7] M. Randeria, in Base-Einstein Condensation, edited by 
A. Griffin, D. Snorke, and S. Stringari, (Cambridge Uni- 
versity Press, Cambridge, 1995). 

[8] J. R. Engelbrecht, M. Randeria, and C. A. R. Sa de Melo, 
Phys. Rev. B 55, 15153 (1997). 

[9] A. S. Alexandrov and S. G. Rubin, Phys. Rev. B 47, 5141 
(1993). 

[10] H. T. C. Stoof, Phys. Rev. B 47, 7979 (1993). 

[11] S. De Palo, C. Castellani, C. Di Castro, and B. K. 

Chakraverty, Phys. Rev. B 60, 564 (1999). 
[12] H. Hu, X.-J. Liu, and P. D. Drummond, Europhys. Lett. 

74, 574 (2006); H. Hu, P. Drummond, and X.-J. Liu, 

Nature Physics 3, 469 (2007). 
[13] M. Y. VeiUette, D. E. Sheehy, and L. Radzihovsky, Phys. 

Rev. A 75, 043614 (2007). 
[14] P. Nikolic and S. Sachdev, Phys. Rev. A 75, 033608 

(2007). 

[15] R. B. Diener, R. Sensarma, and M. Randeria, 
larXiv: 0709.2653 

[16] R. Haussmann, W. Rantner, S. Cerrito, and W. Zwerger, 
Phys. Rev. A 75, 023610 (2007). 

[17] S. Diehl, H. Gies, J. M. Pawlowski and C. Wetterich, 
Phys. Rev. A 76, 021602(R) (2007); jcond^t/0703366j 

[18] See, for example, J. R. Schrieffer, Theory of Supercon- 
ductivity, (Benjamin, Reading, revised printing 1983). 

[19] C. M. Varma, P. B. Littlewood, S. Schmitt-Rink, E. 
Abrahams, and A. E. Ruckenstein, Phys. Rev. Lett. 63, 
1996 (1989). 

[20] We denote by A" = (fc, iuj) and P = (p, iui) collective la- 
bels for momenta and Matsubara frequencies. Integration 
symbols are = {PV)-^ E^,. and = {pvy' E^,^ 
where V is the volume and /3 is the inverse temperature. 
The corresponding normalization of the delta-symbols is 

5k,k' = PVS^^^iS^^k' ■ 

[21] C. Castellani, C. Di Castro, F. Pistolesi, and G. C. Stri- 
nati, Phys. Rev. Lett. 78, 1612 (1997); F. Pistolesi, C. 
Castellani, C. Di Castro, and G. C. Strinati, Phys. Rev. 
B 69, 024513 (2004). 

[22] A. Kreisel, N. Hasselmann, and P. Kopietz, Phys. Rev. 
Lett. 98, 067203 (2007). 

[23] The fields x and (j> as real as functions of space and imag- 
inary time such that in Fourier space X-p ~ Xp and 

4>-p — 4>*p- 

[24] F. Sauh and P. Kopietz, Phys. Rev. B 74, 193106 (2006). 
[25] M. Marini, F. Pistolesi, and G. C. Strinati, Eur. Phys. J. 
B 1, 151 (1998). 



[1] For a recent review see I. Bloch, J. Dalibard, and W. 
Zwerger, larXiv:0704.30TT] 



